Local driving and global interactions in the progression of seizure dynamics 



in 
o 
o 

(N 
5h . 

< 

OO 
(N 

S 1 

Oh' 



CO 

O 

>V 

Oh: 



> 

(N 

(N : 
^t- ■ 
o ■ 
in : 
o . 

o : 
• i— i 

43 : 

Ph. 



X 



Benjamin H. Singer, Miron Derchansky, Peter L. Carlen, and Michal Zochowski 

Neuroscience Program, 
Department of Physics and Biophysics Research Division 
University of Michigan 
Ann Arbor, MI 48109 and 
Division of Cellular and Molecular Biology 
Toronto Western Research Institute, Ontario 
Toronto, Canada 
(Dated: February 2, 2008) 

The dynamics underlying epileptic seizures are well understood. We present a novel analysis of 
seizure-like events (SLEs) in an ex vivo whole hippocampus, as well as a modeling study that sheds 
light on the underlying network dynamics. We show that every SLE can be divided into two phases. 
During the first, SLE dynamics are driven by the intra-network interaction of a network exhibiting 
high internal synchrony. The second phase is characterized by lead switching, with the leading 
region exhibiting low internal synchrony. We show that the second phase dynamics are driven by 
inter-network feedback among multiple regions of the hippocampus. 

PACS numbers: 87.18.Hf, 05.45.Xt, 05.65.Tp 



Epilepsy is a disease that occurs in 2% of the popula- 
tion and is though to be caused by interplay of variety of 
molecular, cellular and network mechanisms lj. However, 
the common underlying feature of all epilepsies is reoc- 
curring seizures due to spontaneous indiscriminate syn- 
chronization and bursting of cell populations 0- Thus, 
it is crucial to understand the fundamental properties of 
seizure dynamics in order to reduce their occurrence. 

The hippocampus is thought to play a central role 
in the genesis of seizures in the most common form of 
epilepsy . While many studies have elucidated the dy- 
namics of seizure onset 0, 0, El j ^ ess m known about 
internal seizure dynamics. Here, we use a novel tech- 
nique to analyze seizure-like events (SLEs) in an ex vivo 
preparation of the whole hippocampus. We have found 
that the observed seizure-like activity can be divided 
into two phases. During the first phase SLE dynam- 
ics are driven by the region of the hippocampus which 
exhibits the highest internal synchrony. Using a com- 
putational model, we propose that in the first phase 
the intra-network dynamics of a single region drive the 
SLE throughout the hippocampus. The second phase is 
marked by the onset of lead switching among regions of 
the hippocampus and, contrary to the first phase, the 
leading region has the lowest internal synchrony. We at- 
tribute this inversion of the relationship between leading 
region and local synchrony to SLE dynamics reflecting 
inter- network feedback among hippocampal regions. 

Seizure activity in animal models of epilepsy has long 
been studied in transverse hippocampal slices 0. How- 
ever, slicing the hippocampus creates artificial two di- 
mensional neuronal circuits which do not capture dy- 
namics resulting from the three dimensional organization 
of the hippocampus. Here, spontaneous seizure-like ac- 
tivity is studied by field potential recording from whole 
hippocampi dissociated from the brains of C57/BL mice, 



as described in The excitability of the neural tis- 
sue is increased by superfusing with low Mg 2+ artificial 
cerebrospinal fluid, thus inducing spontaneous, recurrent 
SLEs. Four electrodes are placed along the temporal- 
septal axis of the hippocampus (Fig. Recordings 
from the four electrodes thus represent neural activity 
in widely separated regions along the longitudinal hip- 
pocampal axis, and reflect both network dynamics within 
the local networks of the four regions, and inter-network 
interactions among them. Bursts recorded in these four 
regions are highly coincident with relatively small lead 
and width variations. Our goal is to characterize how the 
interaction among distinct regions along the temporal- 
septal axis relate to the evolution and termination of in- 
dividual SLEs. 

Recordings from the ex vivo hippocampus reveal 
abrupt transitions from interictal periods to SLEs, which 
are composed of slow bursts and superimposed fast spikes 
(Fig. Q3). To discern the dynamic relationship among 
the four recorded regions, we examine the small tempo- 
ral variations in the onset of slow bursts (Fig.^3, inset), 
as recorded in the four electrodes. We use a modified ver- 
sion of a measure we have ealier developed to quantify 
differences in relative inter-burst intervals (IBIs) between 
electrode pairs (Fig. EJ) , with the underlying idea that ac- 
tivity in a driven region will closely and systematically 
follow the activity of the driving region with minimally 
variable IBIs . The IBIs are calculated separately for 
every electrode pair (i,j) in the the network, and distri- 
butions of IB I values are constructed. The running distri- 
butions of relative IBIs (IBIfj) are updated and renor- 
malized over time with each the onset of each burst. Af- 
ter each update, the Shannon entropy, S = ~^2 k pk^npk, 
of the renormalized IBI distribution is calculated. This 
quantity depends on the relative timing of bursts in both 
electrodes in the pair, and is referred to as a conditional 
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FIG. 1: A: Schematic of the hippocampus showing the temporal 
and septal poles and the approximate placement of electrodes. B: 
Representative local field potential recording of a single SLE (elec- 
trode 2). Inset: a single burst from mid-SLE, showing electrode 
2 (black) leading electrode 1 (gray). C: Two dimensional plot of 
SEj i+SEi j vs. SEj i—SEi j for data shown in B (black). Values 
of SEi j are paired with their nearest temporal neighbors in SEj j. 
Thresholds for SE sum and difference are set based on the distribu- 
tions for a randomly shuffled surrogate dataset (gray). Thresholds 
are set at the 95 th percentile of SEj t i + SEij (horizontal line), and 
the 2.5 th and 97.5 t,l percentiles of SEj^-SE^j (vertical lines). The 
intersection of these thresholds form six sectors. In sector I, activ- 
ity in i significantly leads activity in j, while in sector III, activity 
in j significantly leads activity in i. Events in other sectors are not 
signficant. D: Raster plots showing times of events pairs among ad- 
jacent electrodes {(1, 2), (2, 3), (3, 4)}with significant SE sum and 
difference values. Color denotes septal electrode leading temporal 
(gray) and temporal electrode leading septal (black). The overall 
pattern of leading changes from septal — > temporal (a) to temporal 
— > septal (6) to center — » poles (c). 



entropy of electrode i vs. electrode j (CEij). 

Since the relative IBIs are measured unidirectionally 
(Fig. |2l , the pairwise comparison of CEij and CEji al- 
lows the asymmetric measurement of temporal interde- 
pendencies between the activity at two electrodes. To 
account for the possibility that the temporal interdepen- 
dence between two signals could be an artifact of au- 
tonomous signal properties on one (or both) electrodes 
we calculate the quantity: 



SE hj - 



min(S hi , Sjj) - CE hj 
min(Sia, Sjj) 



(1) 



where Sij is the entropy of (continuously updated) IBI 
distributions observed on a single electrode (IBIi : i) . The 
values of SEs around and below zero indicate that there 
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FIG. 2: The interval IBIij of the i-th electrode with respect to 
the j'-th electrode is calculated as the time difference between an 
event in j and the immediately preceding event in i. CEij is up- 
dated after every calculation of the running distribution of IBIij, 
and assigned to the time of the originating event in i. Analogous 

and Si i. 



is no significant interdependence between the signals, 
whereas if SE tends to one this indicates strong temporal 
interdependence. This measure detects and categorizes 
three basic regimes: the two electrodes are independent 
(both SEij and SEj^ are low); the signal at electrode i 
leads electrode j {SEij is high while SEjj is low); the 
signal is nearly periodic or is completely synchronized 
(both SEij and SEjj are high). 

SE values are computed for every burst of an SLE. 
In order to determine which SE values are due to un- 
derlying neural activity and which are due to chance, 
we compute surrogate datasets by randomly reassigning 
burst times from an individual recording to new elec- 
trode labels (shuffled dataset). Thresholds for signifi- 
cance are then based on the distribution of (SEj^+SEij) 
and (SEj ti - SE itj ) over 10 shuffled datasets (Pig. Hp). 
Events that fall outside the significance threshold for 
both SE sum and difference for each of the electrode pairs 
= {(1,2), (2,3), (3,4)} are depicted in a raster plot 
of individual events and lead pattern (Fig. DP). Within 
the SLE depicted, the focus of activity is initially stable 
and located in the septal (electrode 1) region of the hip- 
pocampus (Fig. Up, region a). The lead pattern rapidly 
switches so that temporal (electrode 4) activity leads sep- 
tal activity (Fig. HP), region b). The lead pattern then 
switches so that a region in the mid-hippocampus (elec- 
trode 2) leads both septal and temporal poles, and this 
pattern remains stable until the SLE ends (Fig. ^3, re- 
gion c). We find that this pattern of septal, temporal, 
and central lead transitions is remarkably stable over dif- 
ferent SLEs recorded from different hippocampi . This 
pattern is observed in 19/23 (82 %) of SLEs analyzed, 
recorded from 6 hippocampi. 

In order to further elucidate this observed switching 
behavior, we build a simple computational model cap- 
turing basic properties of the experimental system. We 
assume that each region recorded by a single electrode 
constitutes a local network which is interconnected with 
other networks along the temporal-septal axis of the hip- 
pocampus. To model the simplest case of two such inter- 
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connected networks, we create two networks of integrate- 
and-fire neurons with a Small World Network (SWN) ar- 
chitecture (N = 15x15 grid each, Fig. EK). The SWN 
architecture constitutes an intermediate type of connec- 
tivity between local and global connectivity [10| , has been 
reported in neural structures, and has been linked to 
seizure generation |(|. The neurons in both networks 
are positioned on a 2-dim lattice with a lattice constant 
a = 1 and periodic boundary conditions. Initially, in each 
network all the neurons within a radius k = 2 are con- 
nected. A fraction (p = 0.3) of those initial connections 
are then randomly rewired. This results in a networks 
having 0.05 connectivity ratio. Additionally, a fraction, 
(/ = 0.3), of randomly chosen neurons in one network re- 
ceives the synaptic current from randomly chosen group 
of Tfi — 10 neurons from the other network. The dynam- 
ics of each neuron are given by: 

^ = - ai v(t) +aJ2 Ja Jik (t) + m (2) 

jec fee/ 

where A = 4 determines the intra-network signal am- 
plitude, B = 0.4 is the inter-network signal amplitude, 
Q-i G (1.0, 1.5) is the membrane leakage coefficient (differ- 
ent for every neuron in the network), and £ £ (0.0, 1.4) 
is a random variable simulating white noise. C denotes 
the set of all neurons connected to i-th neuron via intra- 
network connections, while / denotes the neurons con- 
nected via mter-network connections. is the term de- 
scribing synaptic current arriving from the j-th neuron 
and is given by: 

Jij = (cxp (- — ) - exp {-—)) (3) 

T S Tf 

where t s is the time from the last spike generated at j- 
th neuron; r s = 0.3ms is a slow time constant, whereas 
tauf — 0.03ms is a fast time constant. The interplay of 
those two constants defines the time course of the spike 
decay. 

When the threshold V S pike = 1 is reached, a spike is 
generated and the membrane potential is reset to 0. Dur- 
ing a built-in post-spike refractory period, T = 10ms, the 
membrane does not potentiate in response to incoming 
stimuli. Every neuron in both networks has an additional 
inhibitory mechanism that resets the incoming synaptic 
current to zero if it is below a threshold level (r cu t = 0.3). 
Inclusion of this threshold imposes a requirement for co- 
incident input in spike generation, much as in dendritic 
processing. The parameters if the model (p, A, B, £) are 
set so that both networks are just below a spontaneous 
bursting regime. At a set point in time (t = 10s, Fig. 
|3j3, bottom) T cut is set to zero in network 1, shifting net- 
work 1 into spontaneous bursting. The change in T cut is 
a phenomcnological model of a transiently lowered firing 
threshold, which could be due to multiple neurobiological 
mechanisms. For example, on the single-cell level it could 



be driven by changes in Ca 2+ dynamics and loss of hy- 
perpolarizing activity from Ca 2+ -dependent K + channels 

Bursting in network 1 (Nl) is generated when sufficient 
numbers of neurons spike simultaneously, generating a 
cascade effect in the network. Thus bursting in network 
1 is initially generated through intra-network dynamics 
following a delay after lowering T cut . The activity in 
Nl, in turn, provides input to network 2 (N2), result- 
ing in seizure-like activity in both networks. Bursting 
in N2 is therefore a result of synchronous, inter-network 
signalling, the pattern of which is dictated by the inter- 
nal dynamics of Nl. In this phase of the model seizure, 
when T cut is low, the bursts of Nl lead those of N2 (Fig. 
133, black bar). Nl bursts are also generally higher and 
narrower than N2 bursts, indicating greater internal syn- 
chrony of Nl. 

At t = 60s T cut in Nl is reset to its original value. As 
a result, the dynamics of the 2 network system changes 
dramatically. Bursting dynamics are no longer due to an 
increased intrinsic firing rate of neurons in Nl, but are 
sustained by synaptic input alone. Thus, the bursting 
activity of both networks is not predominantly mediated 
by the internal dynamics of Nl, but by inter-network 
feedback. In this phase, switching in temporal leading 
and internal synchrony are observed among the two net- 
works, due to an inverse driving phenomenon (Fig. 
gray bar). The instantaneous rate of bursting in each 
network is inversely linked to the size and coherence of 
the bursts generated by the other, consistent with burst- 
ing sustained by inter-network interactions For example, 
if a low-amplitude, diffuse burst in Nl, will result in less 
activation of N2. This, in turn, leads to decreased instan- 
taneous frequency and gradual desynchronization of N2. 
Due to the slowing of N2 bursts, Nl activity will tem- 
porally lead bursting in N2. Thus, the less synchronous 
and smaller bursts of one network will lead the more syn- 
chronous and larger bursts of the other until bursting 
ceases. 

To quantify this inverse driving effect, we compare the 
ratio 7 — heig ,^ of bursts in each network to the iden- 

' width 

tity of the leading network for 10 simulated seizures of 
the type shown in (Fig.|23). To measure relationship be- 
tween 7 and temporal leading, we define the expectivity 
of a burst as +1 if the leading burst has a greater value 
of 7 than the following burst, and —1 if the leading burst 
is smaller in 7 than the following burst. In order to cap- 
ture the temporal pattern in expectivity and control for 
variations in burst freqency, we average expectivity over 
a sliding window. In the simulated seizures, expectivity 
is positive while T cut = in Nl, and undergoes a rapid 
transition to negative values when T cut is restored (Fig. 
EK). This finding is consistent with the transition from 
driving by intra-network dynamics to bursting sustained 
by inter-network feedback. 

To determine if similar dynamics may be observed in 
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FIG. 3: A: Schematic connectivity of two connected SWN net- 
works. B: Total synaptic activity from each network. In the first 
phase of the model seizure, activity in Nl leads activity in N2 
(black bar). In the second phase, switching occurs and activity 
in N2 leads activity in Nl (gray bar). Step function denotes the 
period (t G [10,60]) when T cut = 0. 



recordings from the whole hippocampus, we define ex- 
pectivity among two electrodes as: 



E-, 



SE, 
SE, 



> 
< 



(4) 



where W is the burst width. Expectivity (Eij) is calcu- 
lated on a burst-by-burst basis and averaged over sliding 
windows for each adjacent electrode pair. We use burst 
width as a measure of intra-network synchrony, since it 
is less confounded by current source location than burst 
amplitude. In order to account for the variable length of 
SLEs across recordings, we linearly map the time from 
SLE onset to switching onset to the interval [0, 0.5], and 
the time from switching onset to the end of the SLE to the 
interval [0.5, 1.0]. These normalized time series are then 
averaged over 45 comparisons of activity in adjacent elec- 
trodes to calculate average expectivity (E). We observe 
that before the onset of switching (t < 0.5), E > 0, and 
the electrode displaying narrower bursts leads the pair 
(Fig. |3(3, left). After the onset of switching (t > 0.5), 
E < 0, with wider bursts leading (Fig. 0)3, right), as ob- 
served in our model. Thus, these findings are consistent 
with a transition from the intra-network dynamics of a 
single region driving activity throughout the hippocam- 
pus to bursting sustained by inter-network feedback, as 
observed in our model. 
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FIG. 4: A: Average of expectivity among 10 simulated seizures, as 
defined in the text. Left of dashed vertical line: T c ut = 0.0. Right 
of dashed vertical line: Y C ut = 0.3. B: Average (black) ± SEM 
(gray) of expectivity (E) among 45 electrode pairs, as defined in the 
text. Time from SLE onset to switching is normalized to [0,0.5], 
and time from switching to SLE end is normalized to [0.5,1.0]. 
Vertical line: onset of switching. 



To our best knowledge this is the first such record- 
ing, characterization, and proposed mechanism for the 
dynamics of SLE in the whole hippocampus. Using a 
novel measure, we observe and characterize lead switch- 
ing among regions of the hippocampus over the course 
of an SLE. We contruct a simple model to provide in- 
sight into the dynamics underlying these observations. 
We show that each SLE is divided into two phases. The 
first phase is driven by local intra-network dynamics of 
the region with greatest local synchrony, while the sec- 
ond phase is characterized inter-network feedback and an 
inverse leading relationship. 
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